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We study the penguin-dominated B — * PV decays, with P (V) representing a pseudo-scalar 
(vector) meson, in the next-to-leading-order (NLO) perturbative QCD (PQCD) formalism, con- 
centrating on the B — > K<f), nK* , pK, and cuK modes. It is found that the NLO correc- 
tions dramatically enhance the B — » pK, uK branching ratios, which were estimated to be 
small under the naive factorization assumption. The patterns of the direct CP asymmetries 
A C p(B° -> p^K*) « A C p(B ± -> p K±) and \A C p(B° -> ir^K**)] > \A C p(B ± -> tt J^* ± )| 
are predicted, differing from \A C p(B° -> tt t /<"*) > \A C p(B ± -» tt ^)]. The above patterns, if 
confirmed by data, will support the source of strong phases from the scalar penguin annihilation in 
PQCD. The results for the mixing-induced CP asymmetries Sf are consistent with those obtained 
in the literature, except that our S p o K is as low as 0.5. 
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I. INTRODUCTION 

It has been a challenge to explain the branching ratios of the penguin-dominated B — > PV decays, such 
as B — > K(f>, wK* , pK, and ujK. The predictions from the naive factorization assumption (FA) [lj and from 
the QCD-improved factorization (QCDF) approach 0] are usually smaller than the observed values [!, 0, [5[ . 
In order to account for the data in QCDF, the hadronic parameters associated with the penguin annihilation 
amplitudes must be tuned (i.e., adopting the scenario "S4"), so that these amplitudes are large and constructive 
to the leading FA ones @ . On the contrary, the B — > K<j> branching ratios have been predicted correctly @, Q in 
the perturbative QCD (PQCD) approach [EG3- The reason is that the leading amplitudes are not factorizable 
in QCDF based on collinear factorization theorem due to the existence of end-point singularities. Hence, they 
are characterized by a scale of order mj, ms being the B meson mass. The characteristic scale of the leading 
amplitudes in PQCD, being factorizable since they are formulated in kx factorization theorem, is of order 
•v/mflA, where A is a hadronic scale. At this lower characteristic scale, the Wilson coefficients of the QCD 
penguin operators are enhanced, so that the PQCD predictions for the B — > K<j> branching ratios are large 
enough. This mechanism has been called as the dynamical enhancement of penguin contributions @ . 

The PQCD analysis of the B — > K<f> decays in 0, Q was performed at leading order (LO). Recently, the 
PQCD formalism has been extended to next-to-leading-order (NLO) accuracy where the contributions 
from the NLO running of the Wilson coefficients, the vertex corrections, the quark loops, and the magnetic 
penguin were taken into account. It was found that these NLO contributions affect penguin-dominated decays 
more than tree-dominated ones. For example, the quark-loop and magnetic-penguin corrections decrease the 
penguin amplitudes by about 10%, and the B — > ttK branching ratios by about 20% However, these NLO 
contributions have almost no impact on the tree-dominated B — > irir branching ratios. Therefore, it is essential 
to examine the NLO effects on LO predictions for other penguin-dominated modes. In this paper we shall study 
the penguin-dominated B — > PV decays in NLO PQCD, concentrating on B — > K<fi, irK*, pK, and ujK . It is 
expected that the LO PQCD predictions for the tree-dominated modes, such as B — > pir, urn [l2j], are stable 
under the above NLO corrections. 

The second motivation to study the B — > PV decays arises from the observed direct CP asymmetries of the 
B ± — > p^K^ decays: the data of Aqp{B ± — * 7r°AT ± ) almost vanish, but those of Acp{B ± — > p°K ± ) are sizable 
[T3 |. The direct CP asymmetries of both the charged B meson decays were predicted to be large in LO PQCD 
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However, as demonstrated in [Tl|, the aforementioned NLO corrections can lower the LO predictions for 
\Acp{B ± — + tt° K^)\ to the measured values. The responsible mechanism comes from the vertex corrections, 
which induce a large and almost imaginary color-suppressed tree amplitude. This amplitude then decreases 
the relative strong phase between the total tree amplitude and the total penguin amplitude involved in the 
_B ± — > ifiK^ modes. Hence, it is natural to investigate whether the LO predictions for Acp(B ± — > p°K ± ) and 
for the direct CP asymmetry of another charged mode, Acp{B ± - ► ^K*^), will be modified by the vertex 
corrections significantly. 

The third motivation concerns the extraction of the standard- model parameter 8111(20!), cf>i being the weak 
phase, from the mixing-induced CP asymmetries Sf of the penguin-dominated B — > / decays. In principle, 
the results should be identical to those from the tree-dominated b — > ccs transitions, such as B — > J/tpK^*\ 
Nevertheless, potentially substantial deviations AS = Snepgaip — Sr.™ < between the above two extractions 
have been observed [13]. Theoretical estimates [HI EE [iEll3 til, H S, HH, H3, HI] of the tree pollution in the 
penguin-dominated decays arc then crucial for justifying these discrepancies as promising new physics signals, 
which, however, gave tiny and positive AS for most modes so far. The mixing-induced CP asymmetry of the 
B° — > ir°Ks decay has been calculated in NLO PQCD [ll|], and the result is consistent with those obtained in 
the literature. It is necessary to complete the survey by computing the mixing-induced CP asymmetries of the 
penguin-dominated B — > PV decays. Since the deviations in the penguin-dominated modes are mainly caused 
by the color-suppressed tree amplitude, which is sensitive to the vertex corrections, we shall adopt the NLO 
PQCD formalism here. 

It will be demonstrated that the dynamical enhancement of the penguin amplitudes in the B PV decays 
sustains the NLO corrections. Especially, the NLO effects dramatically enhance the B — » pK, ujK branching 
ratios, which were predicted to be small in LO PQCD 14] and in FA @|. The NLO PQCD results for all 
the considered B — > PV branching ratios then match the data within theoretical uncertainties. Note that the 
measured B — > luK branching ratios are close to the B — » pK ones, though the former receive the additional 
contribution from the Wilson coefficient 05. Because the magnitude of 05 increases rapidly at a low energy scale, 
the similarity between these two modes implies that two-body nonlcptonic B meson decays are insensitive to 
low-energy dynamics. Their data, if persisting, will encourage the application of a perturbation theory to these 
processes. 

The LO PQCD prediction for the direct CP asymmetry Acp{B ± — ► p°K ± ) [l4| turns out to be stable under 
the NLO corrections. That is, we predict the different patterns Acp(B° — > p T K^) w AcpiB^ — > p°K ± ) and 
\Acp(B° — > 7r =F A' ± )| 3> \Acp{B ± — > 7r°A ± )|, which are attributed to the different orientations of the complex 
penguin amplitudes: the penguin amplitude is inclined to the imaginary (real) axis in the former (latter). When 
the penguin amplitude is almost imaginary, its strong phase relative to the tree amplitude is maximal. The 
color-suppressed tree amplitude, despite of being enhanced by the vertex corrections, does not change it much. 
The configuration of the relevant amplitudes in the B — > 7r K* decays is somewhat located between the B — ► pK 
and B -> irK ones, giving \A C p(B° -> tt^ K*^)\ > \A C p(B ± -> 7r if* ± )|. We point out that the above patterns 
can not be attained by the other theoretical approaches, including QCDF soft-collinear effective theory 

(SCET) [H[2l|,[|l, and QCDF plus final-state interactions (FSI) [H[. 

As to the mixing-induced CP asymmetries, we obtain AS* > at LO for all the B — > PV modes. The NLO 
corrections increase most of AS except that for B° — ► p°K$, whose mixing-induced CP asymmetry becomes as 
low as 0.5. The uniqueness of the B° — ► p°Ks decay results from the fact that both the penguin amplitude and 
the color-suppressed tree amplitude are almost imaginary and parallel to each other at NLO, a configuration 
which leads to the very negative AS*. Among the penguin-dominated decays we have investigated (B° — > ir°Ks, 
4>Ks, p°Ks, uiKs), B° — * 4>Ks is the cleanest one with the minimal tree pollution for extracting sin(2</>i). 
Roughly speaking, the NLO PQCD predictions for the B — > PV mixing-induced CP asymmetries are consistent 
with those in QCDF [13, El, but differ from those in QCDF plus FSI [l|, which gave AS w +0.04 for the 
B" -> p°K s decay 0. 

We review the LO PQCD factorization formulas for the B — > PV decays, and derive the NLO ones in Sec. II. 
The branching ratios, the direct CP asymmetries, and the mixing-induced CP asymmetries of the penguin- 
dominated B — > PV modes are calculated using the NLO PQCD formalism in Sec. III. We have updated the 
meson distribution amplitudes and the relevant standard-model inputs in the calculation. Section IV is the 
conclusion. 
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II. FACTORIZATION FORMULAS 

The amplitude for a B meson decay into the two-body final state M 2 M 3 through the b — > s transition has 
the general expression, 



a(b - m 2 m 3 ) = v; b v us + v; b v cs a[% M3 + v* b v ts A% M3 
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where fs (/m 2 , Jm 3 ) is the B meson (M 2 meson, M 3 meson) decay constant, F e (M e ) the color-allowed fac- 
torizable (nonfactorizable) tree emission contribution, F c m 3 (M £ m 3 ) the color-suppressed factorizable (nonfac- 
torizable) tree emission contribution, F a (M a ) the factorizable (nonfactorizable) tree annihilation contribution, 
and those with the additional superscripts P the contributions from the penguin operators. With the amplitude 
A(B — ► M 2 M 3 ) in Eq. {l}, we derive the branching ratios, the direct CP asymmetries, and the mixing-induced 
CP asymmetries of two-body nonleptonic B meson decays, which have been defined in [ll[. 
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TABLE I: B -> K<f> decay amplitudes in LO PQCD. 



The LO PQCD factorization formulas for the B — > i4T</), ^Jf*, and p(co)K decay amplitudes are summarized 
in Tables HI HH and lllli respectively. The Wilson coefficients for the factorizable contributions and a'^ for 
the nonfactorizable contributions are the combinations of in the standard definitions [xij . where q = u, d or 
s denotes the quark pair produced in the electroweak penguin. The explicit expression for each amplitude in 
Tables U tillll is similar to that of the B — > P 2 P 3 modes [1JJ], but with the replacements, 

0^ 03 , <f>3 ~> 03 > ^ ^ 4 . m 02 -* "7,02 , m 03 -> -m 3 , (3) 
for B — > P2V3 (involving the B — > P 2 transition), and 

02 -> 02 , 0^ 02 1 02 02 ) m 02 -m 2 , -> -m 03 , (4) 

for B — ► V 2 P 3 (involving the P — > V 2 transition), where mo is the chiral enhancement scale associated with a 
pseudo-scalar meson, and m the vector meson mass. The above replacements also apply to the factorization 
formulas for the NLO corrections from the quark loops and from the magnetic penguin. 

We add the various NLO corrections to the LO factorization formulas in Tables Q] IIII1 The NLO corrections 
to the smaller nonfactorizable contributions will be neglected. It is understood that the Wilson coefficients 
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appearing below refer to the NLO ones. The vertex corrections to the B — ► PV decays modify the Wilson 
coefficients for the emission amplitudes in the standard definition into 

«»00 - a 2 (M) + ^ F ^y 2 (M), 

a,( M ) - + ^C F ^^V t (M) , i = 3 - 10 , (5) 

where M denotes the meson emitted from the weak vertex, and the upper (lower) sign applies for odd (even) i. 
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TABLE II: B -> ttF"* decay amplitudes in LO PQCD. 
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When M is a pseudo-scalar meson, the functions Vi(M) are given, in the NDR scheme, by @, Q 



12 In— -18 + -^ — - 



IM Jo 



dx(f>f I (x)g(x) , 



for i = 1 - 4, 9, 10 



, to 6 2a/2A^ 
-12 In— + 6 



jti 

Im Jo 



Sm 



dx 4>m{x) g(l — x) , for i = 5, 7 



4>m(x) h(x) 



for i = 6, 8 



(6) 



a; being a momentum fraction. For a vector meson M, (fi^ (</>£f) is replaced by 4>m [~4 >s m)i an d /m by /^in 
the third line of the above formulas. The explicit expressions of the hard kernels g and h can be found in [6(. 
When a vector meson is emitted from the weak vertex in, for example, the cases of B — > if <fi, irK*, the leading 
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TABLE III: B — > pF, ojF decay amplitudes in LO PQCD with the upper (lower) signs applying to the pK (u)K) modes. 
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vertex corrections, the NLO piece Ogy C , containing the second terms of ae,8 in Eq. contributes through the 
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These additional amplitudes have emerged in the NLO analysis of the B — > V V decays [13] . 

Taking into account the NLO contributions from the quark loops and from the magnetic penguin, the decay 
amplitudes in Eq. ([I]) are modified into 
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where the upper (lower) signs apply to the pK (loK) decays. The NLO amplitudes M.^ M ' -^M Af ' -^M M ' 

and A^j^ Ms denote the up-loop, charm- loop, QCD-penguin-loop, and magnetic-penguin corrections, respec- 
tively. Their expressions are similar to those for the B — > PP decays [TTj . but with the replacements specified 
in Eqs. ([3]) and ((4]). The magnetic-penguin contribution to the B — > PV modes has been computed in (28|. It 
will be shown that the quark-loop corrections are always constructive, increasing all the B — > PV branching 
ratios considered here. The magnetic-penguin corrections decrease the B — > K<p, ttK*, ojK branching ratios, 
but enhance the B — ► pK ones slightly. 



III. NUMERICAL RESULTS 



We perform the numerical analysis in this section based on the factorization formulas derived above. The 
choices of the B meson wave function, the B meson lifetimes, the decay constants and the chiral enhancement 



scales associated with the pseudo-scale mesons, and the weak phase 4>\ — 21.7° are the same as in We take 
Trip — 0.77 GeV, m u = 0.78 GeV mx* = 0.89 GeV, and = 1.02 GeV for the vector meson masses (29|, and 
rrib = 4.8 GeV for the b quark mass appearing in the magnetic-penguin operator 0% g . The longitudinal decay 
constants of the vector mesons can be extracted from other decay rates, for example, f p ,K* from the measured 
t~ — » (p~ , K*~)v T decay rate s l29l.l30ll. The transverse decay constants have been evaluated in QCD sum rules, 
which are summarized below |30l. l3ll. l32j: 

f p = (0.209 ± 0.002) GeV , /J = (0.165 ± 0.009) GeV , 

fu = (0.195 ± 0.003) GeV , /J = (0.145 ± 0.010) GeV , 

f K , = (0.217 ± 0.005) GeV , fj { , = (0.185 ± 0.010) GeV , 

U = (0.231 ± 0.004) GeV , /J = (0.200 ± 0.010) GeV . (9) 
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We employ the updated twist-2 and twist-3 pseudo-scalar meson distribution amplitudes, and the updated 
twist- 2 vector meson distribution amplitudes from QCD sum rules [3l|, [HI]. For the twist-3 vector meson 
distribution amplitudes, we adopt the asymptotic models [34] ]. since some of them are still uncertain in sum 
rules. The explicit expressions of the above distribution amplitudes are collected in the Appendix. The sum- 
rule analysis of the kaon and K* meson distribution amplitudes have been also performed in [35| and [36[, 
respectively. When analyzing the theoretical uncertainty, we consider the following ranges for the Gegenbauer 
coefficients in the light-meson distribution amplitudes, whose definitions are referred to the Appendix, 

al = 0.25 ±0.25 , af = 0.25 ± 0.25 , 

a\ p = al u = 0.15 ±0.15 , al K , = 0.11 ± 0.11 , a!j = 0.0 ± 0.2 . (10) 

Since our results do not vary with the Gegenbauer coefficients af and a\ K , much, their values have been fixed at 

af = 0.06 and a{ K , = 0.03. All the Gegenbauer coefficients associated with the twist-3 distribution amplitudes 
are also fixed for simplicity. We update the CKM matrix elements |29J , 

V ud = 0.97377, V us = 0.2257, \V ub \ = (4.31 ± 0.30) x HT 3 , 
V cd = -0.230 , V cs = 0.957 , V cb = 0.0416 , 

03 = (70 ±30)°. (11) 

The increase of |T4f>| from its old value (3.67±0.47) x 10 -3 will enhance the direct CP asymmetries of B meson 
decays as shown later. 

We then evaluate the relevant transition form factors at maximal recoil, 

F^ = 0.271832 , = 0.4318;8 , 

= 0.32t°;°£ , A^ = 0.2918:88, A ^=0.38l8:8L (12) 

where the theoretical uncertainties come from the variation of the parameters associated with the B meson 
wave function, and from Eqs. (UHl) and (HH). The above values of F^ and F S K are higher than F B n = 0.24 
and Fq K = 0.36 obtained in [ll| from the old pseudo-scalar meson distribution amplitudes. The others are 
similar to A^ p = 0.30, = 0.28, and A® K " = 0.37 from QCD sum rules [30:]. Our A® K * is a bit above that 
from the covariant light-front QCD calculation [37J. Note that larger values of the B — > V form factors have 
been adopted in QCDF Q in order to account for the observed B — * PV branching ratios. 
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TABLE IV: RG evolution of the Wilson coefficients involved in the B — > uiK decays, where a; (At) and cnvc{n) represent 
the first and second terms on the right-hand side of Eq. ([5]), respectively. 

We present in Table IIVI the renormalization-group (RG) evolution of the Wilson coefficients at and of the 
vertex corrections aivc i n the scale fi for the B — > luK decays, which correspond to the first and second terms 
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on the right-hand side of Eq. respectively. The RG evolution exhibited by the other B — ► PV modes is 
similar. It is seen that the coefficients a4_@ associated with the penguin operators are enhanced more than the 
leading tree coefficient a\ at a low fi. This mechanism, referred to the dynamical enhancement in gives the 
large branching ratios of penguin-dominated B meson decays 0, H El ■ Most of the vertex corrections remain 
subleading even at a low /i, except those to 0,2 and 0,5. Though the correction a^vc to 0,3 is significant, its 
effect, much smaller than that of a^vCi is negligible. The large a^yc has been found to be responsible for the 
reduction of the LO PQCD predictions for the direct CP asymmetry \Acp{B ± — > ^K^)] [ll[ as explained in 
the Introduction. The new effect observed in this work arises from the large a§vc> which will be elucidated 
below. 
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19.8 e~ i0 - 5 


P' 


11.1 e* 


31 11.4e-' 31 


5.8 e 4 


31 5.6 e"* 3 


9.9 e l 


3.1 10 .i e l31 


3.1 e 4 


3.1 


3.2 e* 31 



TABLE V: Topological amplitudes in units of 10" 5 GeV. 

Since all the meson distribution amplitudes derived from QCD sum rules are defined at the scale 1 GeV, we 
propose to freeze the RG evolution of the Wilson coefficients at this scale, when the energy runs below it. That 
is, we set Ci(fi) = Q(l GeV) as \i < 1 GeV. Note that the freezing scale was chosen as 0.5 GeV in [O, We 
have checked that the new choice of the freezing scale decreases the NLO PQCD predictions for the B — > wK 
branching ratios by about 15%, and has a tiny effect on those for the B — * tttt branching ratios. The final 
outcomes are still within the theoretical uncertainty in [ll|, [2tJ ■ In order to facilitate the discussion, we present 
the central values of our predictions in terms of the topological amplitudes in Table [V] where T', C", P', P cl , 
and P' ew denote the color-allowed tree, color-suppressed tree, QCD penguin, color-suppressed QCD penguin, 
and electro- weak penguin contributions, respectively. Their definitions are referred to [3g|. Note that P cl exists 
only in the B — > Kef) and B — > uK decays due to the flavor-singlet quark content of the <f> and u> mesons. We 
have estimated the averaged hard scale (t) for each of the above topological amplitude, and found (t) f» 2.0 
GeV for T and P' ew , (t) » 1.5 - 2.0 G eV for P', and (t) » 1.0 - 1.5 GeV for C and P cl . That is, the hard 
scales, despite of being all of 0(\/tobA), vary a bit among the different topological amplitudes. 

The NLO PQCD results for the B — > PV branching ratios are displayed in Table IVI1 where the QCDF 
predictions from the scenario S4 for the hadronic parameters 0] are shown for comparison. It has been known 
that QCDF gives the branching ratios similar to those in PQCD, if adopting S4. The QCDF predictions from 
the default scenario Q are usually smaller than the data by a factor 2. It is found that the NLO Wilson 
evolution, labelled by LOnlowCj has a more significant effect on the B — * Kcf> decays than on the B — > ttK 
ones 11]. The reason is that the former depends on the Wilson coefficient 05 (see Table HJ, which is enhanced 
more at NLO compared to etg in the latter. However, 05 also receives a significant destructive vertex correction 
as shown in Table IIV1 such that the B — > K<p branching ratios drop rapidly from the column LOnlowc to 
+VC in Table IVIl The magnetic penguin correction, labelled by +MP, further decreases the branching ratios 
by about 30%. Eventually, the B — > K<fi branching ratios in NLO PQCD are still consistent with the measured 
values. 

For the B — > irK * branching ratios, our values in the column LOnlowc are close to the LO PQCD predictions 
in [3^ |. The effects from the vertex corrections, the quark loops, and the magnetic penguin follow the pattern 
appearing in the B — * nK modes flll | . The main difference between these two decays is that the former involve 
only a4 (no ag). Therefore, the B — ► nK* branching ratios are expected to be much smaller than the B — > ttK 
ones. As indicated in Table [VTl our NLO PQCD predictions (central values) for the former are about 1/3 of 
those for the latter. Confronting with the data, the NLO predictions for the B ± — > 7r ± K*°, i? ± — > ti^K^, 
and B° — > 7r T A"* ± branching ratios reach only 2/3 of the measured values. Nevertheless, considering the 
uncertainties of both the theoretical predictions and the experimental data, the discrepancy is not serious. On 
the contrary, the PQCD results for the B° — > ir°K*° branching ratios are in good agreement with the data. 
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Mode Data [13] QCDF(S4) [6] LP LO NLO wc +VC +QL +MP +NLO 



B ± 


-> K ± <f) 


8.30 ±0.65 


11.6 


13.8 


41.8 


14.3 


44.3 


28.4 


7 8+ 5 ' 9 
'•8-1.8 


(+5.6) 
(-1.7) 


B° 


-> K°<t> 


8-31^ 


10.5 


12.9 


39.2 


13.4 


41.5 


26.6 


7 q+5.4 

'■°-1.6 


(+5.1) 
(-1.5) 


B ± - 


-4 ^K'° 


10.7 ±0.8 


8.4 


5.5 


10.6 


9.6 


11.3 


6.5 




(+2.7) 
(-1.4) 


B ± - 




6.9 ±2.3 


6.5 


4.0 


6.8 


6.2 


7.1 


4.5 


4.3+H 


(+1-7) 
(-1.0) 


B°- 


-> n^K** 


9.8 ± 1.1 


8.1 


5.1 


9.3 


8.8 


9.9 


6.0 


6.0+- 


(+2.4) 
(-1.3) 


B° - 


-k°K*° 


1.7 ±0.8 


2.5 


1.5 


3.5 


3.3 


3.8 


2.0 


2 +1 ' 2 


+0.9) 
(-0.4) 


B ± 
B ± 




< 48 


9.7 
4.3 


3.6 
2.5 


6.1 
4.0 


7.8 
5.0 


6.4 
4.3 


6.7 
3.8 




(+6.4) 
(-4.3) 
(+3.6) 
(-2.6) 


B° 




9-9±i:S 


10.1 


4.7 


7.0 


7.8 


7.2 


7.9 


8.8:11 


(+6.2) 
(-3.9) 


B° 




5.6 ± 1.1 


6.2 


2.5 


3.5 


3.9 


3.5 


4.4 




(+3.2) 
(-2.0) 


B ± 
B° 


-> wif* 
-a;* 


6.9 ±0.5 
4.8 ±0.6 


5.9 
4.9 


2.1 
1.9 


5.9 
6.0 


9.3 
8.6 


6.4 
6.6 


4.8 
4.8 


10.6+f 8 ' 

9 8+ 8 ' 6 


t(+7.2) 
(-4.4) 
(+6.7) 
(-4.3) 



TABLE VI: Branching ratios in the NDR scheme in units of 10 . The label LOnlowc means the LO results with the 
NLO Wilson coefficients, and ±VC, ±QL, ±MP, and ±NLO mean the inclusions of the vertex corrections, of the quark 
loops, of the magnetic penguin, and of all the above NLO corrections, respectively. The errors in the parentheses arise 
only from the variation of the hadronic parameters. 

Our LO PQCD results for the B — » pK decays in Table IVT1 are close to those obtained in [l4|]: B(B^ — > 
p ± K°) = 2.96 x 1(T 6 , B{B ± -> p ^) = 2.18 x 10" 6 , B(B° -> p^K^) = 5.42 x 10~ 6 , and B(B° -> p°K°) = 
2.49 x 10~ 6 . That is, the LO predictions fall short by a factor 2 compared to the data [HI]. The B — > pK 
branching ratios, also smaller than the B — > nK ones, are attributed to the destructive combination of the 
Wilson coefficients 04 — 2(m§K I fnB)&% in the former, and to the constructive combination ± 2(moK / 'mfi)a6 
in the latter [see the replacement rriQK — > —n^oK in Eq. Q]. The NLO Wilson evolution enhances the B — ► pK 
branching ratios, but not sufficiently. Since the destructive combination flips sign of the penguin contribution 
to the B — > pK decays (comparing Table [V] here and Table V in [HI]), the NLO effects change accordingly. For 
example, the magnetic-penguin correction becomes constructive. The B — > pK branching ratios then increase 
all the way up to the measured values as shown in Table IVIl It turns out that the relation B(B — > pK) < 
B(B -► ttK *) in LO PQCD and in FA Q is reversed into B(B -> pK) > B(B -> irK*) at NLO. 

Our LO PQCD results for the B — > loK branching ratios are also consistent with those in 1141 : B(B ± — * 
wK*-) = 3.22 x 1(T 6 and B(B° ujK°) = 2.07 x 1(T 6 , which are smaller than the data [lj Similar to 
B — > i^t/), the i? —> wif decays involve the Wilson coefficient 05 through the color-suppressed QCD penguin 
amplitude F^ K , in addition to 04 and (see Table ITTT1) . Hence, the enhancement of the B — > loK branching 
ratios from the NLO Wilson evolution is substantial. Because of the sign flip of the QCD penguin contribution 
as in B — > p°K, the vertex correction to 05 becomes constructive, opposite to that in the B — > K(j) case (see 
the opposite contributions of P cl in these two modes in Table IVf . This explains the jump of the B — > uoK 
branching ratios from the column LOnlowc to ±VC in Tabic IVIl Though the magnetic penguin correction is 
destructive, the net NLO predictions for the B — > loK branching ratios remain large, with the central values 
being higher than the observed ones. 

As stated before, 05 and its associated vertex correction exhibit dramatic running effects at a low scale [lT| . 
Therefore, if lifting the hard characteristic scales in the factorization formulas slightly (notice (t) ~ 1.0—1.5 
GeV for P cl ), the 05 contribution will be moderated. Then the predicted B — > Kcj> and B — > uiK branching 
ratios will increase and decrease, respectively, approaching the central values of the data. However, we shall 
not attempt such a fine tuning here, but point out that the comparison of the two branching ratios provides 
an interesting test on the significance of the as contribution, i.e., on the sensitivity of two-body nonleptonic B 
meson decays to low-energy dynamics. Note that the relevant form factors Aq P and Aq^ are roughly equal as 
indicated in Eq. p2[) . If allowing the hard scale to be as low as 0.5 GeV, B(B — + ujK) would become three 
times of B(B — > p°K). Because the measured branching ratios are close to each other, as plays a minor role, 
and a higher hard scale is implied. It is then likely that two-body nonleptonic B meson decays are insensitive 
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to low-energy dynamics, and a perturbation theory is applicable to these processes. 



Mode Data [13] QCDF(S4) [6] LP LOnlowc +VC +QL +MP +NLO 



B ± -> K^fy 
B° -> K°4> 


3.7 ±5.0 
9 ±14 


0.7 
0.8 


-2 



-1 



-1 




1 


-1 



1+0 (+0) 

-i(-i) 

i'_i_n~i 
o i i In-vj 


B ± -> iv+K* 

B — > 7T K 

B° -> ^iT* 


-8.6 ±5.6 
4 ± 29 
-5 ±14 

1+27 

-26 


0.8 
—6.5 
-12.1 

1.0 


-3 
—38 
-56 
-5 


-2 
—31 
-40 
1 


-2 
—21 
-42 
-13 




—29 
-38 
2 


-2 
—45 
-60 
3 


i+i (+i) 

i -O(-0) 
o 9 +21(+16) 
13 -28 (-19) 
fif) +32 (+20) 
UU -19(-15) 
11 +7 (+5) 
11 -5(-2) 


B± - p±A° 
B± -> ± 
B° -> p T A± 
B° -» p°if° 


Q1+H 

ol -10 
1? +15 
± '-16 


0.8 
31.7 
20.0 
-2.8 


2 

79 
83 
7 


1 

74 
72 

-8 


1 

73 
69 
3 


1 

71 
71 

-8 


1 

79 
64 

-6 


1±1(±1) 
71 +25(+l7) 
' -38 (-14) 
g 4 +24(+7) 
" -30 (-11) 
7 +8 (+7) 
-5 (-4) 


B± -> wJS^ 
B° -> ^" 


5±6 


19.3 
3.7 


82 

-4 


45 
7 


37 

-2 


43 
7 


57 

8 


o 9 +15(+4) 
JZ -17(-5) 
n+2(+2) 
J -4(-3) 



TABLE VII: Direct CP asymmetries in the NDR scheme in percentage. 

The NLO PQCD predictions for the direct CP asymmetries of the penguin-dominated B — > PV decays are 
listed in Table IVlTl where the QCDF results from S4 Q are also shown. The QCDF results from the default 
scenario are always opposite in sign compared with the data. Our predictions for the direct CP asymmetries have 
larger magnitude here, since the updated CKM matrix element \V u b\ has increased substantially. If adopting 
\V u b\ extracted from exclusive decays alone, and a lower B meson decay constant, the magnitude of the above 
direct CP asymmetries could drop by 30% easily. To be cautious, we emphasize only the various patters of the 
direct CP asymmetries exhibited by the considered B — > PV decays, instead of their actual values. Because 
the NLO Wilson evolution increases the penguin amplitudes, i.e., the B — ► PV branching ratios, it dilutes the 
direct CP asymmetries. The effects from the quark-loop and magnetic-penguin corrections can be understood 
in the same way. The direct CP asymmetries of the B — > K<f> modes, which are almost pure-penguin processes, 
remain vanishing at NLO as shown in Table fVTTl 

We explain the different effects from the vertex corrections on the direct CP asymmetries Acp{B^ — ► Tt^K*^), 
A-cpiB^ — ► p°K^), and Acp(B ± — > ir^K^). It has been known [ll[ that the penguin amplitude P' in the 
B — > ttK decays is in the second quadrant, with a strong phase about 15° relative to the negative real axis, 
and that the source of this phase arises mainly from the almost imaginary scalar penguin annihilation. The 
color-allowed tree amplitude T" is roughly aligned with the positive real axis. Consequently, the B° — ► n^K^ 
decays, involving P' and T", show a sizable direct CP asymmetry @. On the other hand, the color-suppressed 
tree amplitude C , enhanced by the vertex corrections, becomes almost imaginary. It then orients the sum 
T 1 + C into the fourth quadrant, such that T" + C and P' + P' ew more or less line up (point to the opposite 
directions) in the B ± — * n°K ± decays. This is the reason the magnitude of Acp{B ± — > -k^K^ tends to vanish 
in NLO PQCD, leading to the pattern [ill, 

\A CP (B ± -> n^K^l » \A CP (B ± -> 7r°i^ ± )| , (13) 

which is in agreement with the data. 

The pattern of the direct CP asymmetries in the B — > irK* decays can be described similarly, but with 
the distinction that their penguin emission amplitudes contain only 04 (no ae). The almost imaginary scalar 
penguin annihilation then renders P' more inclined to the positive imaginary axis compared to that in the 
B — > ttK case. Hence, T 1 and P' are more orthogonal to each other in the B° — > tt^K*^ 1 decays. For this 
configuration, the modification from C will be weaker. That is, the total tree amplitude T' + C and the total 
penguin amplitude P' + P' ew will not line up completely in the B*- tt^K*^ modes. We then predict the larger 
B — * ttK* direct CP asymmetries, and the pattern shown in Table FyTA 

\A CP (B ± - n^K^)] > \A CP (B ± -> Tr ^^)! . (14) 
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The situation in the B — ► pK decays is even more extreme, where the real part of P' almost diminishes 
due to the destructive combination of the Wilson coefficients 04 — 2(uiok /msjog. Because of the sign flip of 
the penguin annihilation amplitude in Table ITTT1 under the replacements — ► — m p and m^K - * ~moK 
in Eq. (j4|), P' is roughly aligned with the negative imaginary axis as indicated in Table Ivl This explains the 
predicted sign of Acp(B° — > p T K^) opposite to that of Acp(B° — > tt^K^). In this case T" and P' have the 
maximal relative strong phase, such that the modification from C is not obvious. We then predict the larger 
B — > pK direct CP asymmetries, and the pattern, 

A CP {B± - p*K±) » Acp{B± -> p ^) . (15) 

That is, the NLO corrections to the P — ► p-RT direct CP asymmetries are minor, and the LO predictions derived 
in [l4| are reliable. 

The direct CP asymmetries in the B — > wif decays can be understood by means of those in B — > pfsT plus 
the effect from the color-suppressed QCD penguin amplitude P cl . Since P cl is roughly aligned with the positive 
real axis as shown in Table [Vl the relative strong phase between the tree and penguin amplitudes decreases 
from the maximal configuration in the B — > pK decays. Therefore, the vertex corrections have some effect on 
AcpiB^ — > ujK^), which is then expected to be smaller than AcpiB^ — > p°K ± ), but still sizable as indicated 
in Table IVHl We stress that the above predictions for the direct CP asymmetries, if confirmed by the future 
data, will support the source of strong phases from the scalar penguin annihilation in PQCD. 



Data [13] LO LO NL owc +VC +QL +MP +NLO 

S^ks 0.47 ± 0.19 0.72 0.72 0.72 

S p o Ks 0.17 ± 0.58 0.71 0.66 0.47 

S uKs 0.64 ± 0.30 0.76 0.69 0.84 



TABLE VIII: Mixing-induced CP asymmetries. 

The NLO PQCD predictions for the mixing-induced CP asymmetries Smk S i M — <fi, p°, and u>, are collected 
in Table IVIIII It is found that all Smk s exhibit positive deviations from those of the tree-dominated b — > ccs 
transitions at LO. The deviation in the B° — ► (pKg decay remains negligible at NLO due to the absence of 
the tree pollution from the color-suppressed tree amplitude. It is interesting to see that the NLO corrections 
turn S p ox s into a much lower value about 0.5. The reason is, again, the destructive combination of the Wilson 
coefficients 04 — 2(mo_fs"/ms)a6, which flips the sign of the penguin amplitude P' as shown in Table [Vl The 
color-suppressed tree amplitude C", modified by the vertex corrections, then becomes parallel to P' , leading 
to the very negative deviation. This observation is consistent with the tendency indicated by the data. The 
deviation of S^Ks 1S positive, since the sign of P' is opposite to that in the B° — > p K$ decay as shown in 
Tables ILTll and IVl The large deviations exhibited by the B° — > p°Ks, loKs decays at NLO are attributed to the 
involved large C (see Table [Vj . Our predictions are basically in agreement with those in QCDF [Tt], EH 1 but 
differ from those in QCDF plus FSI [19j], which predicted S p o Ks 0.76 with a positive deviation. 

IV. CONCLUSION 

In this paper we have surveyed the penguin-dominated B — ► PV decays in the NLO PQCD formalism, 
concentrating on the B — ► K<fi, irK*, pK, and cuK modes. The NLO Wilson evolution enhances all the 
branching ratios, while the other NLO corrections modulate the branching ratio of each mode in a different way. 
For example, the vertex corrections to the Wilson coefficient 05 decrease (increase) the B — > Kef) (B — > loK) 
branching ratios significantly. The net NLO effects modify the LO predictions for the B — > K<f), ttK* branching 
ratios a bit, but increase those for the B — > pK, loK ones by a factor more than 2. Owing to this dramatic 
enhancement, the NLO PQCD results for these B — > PV branching ratios are in better agreement with the 
data. Since the penguin emission amplitudes in the B — > irK, ttK*, and pK decays are proportional to 
04 ± 2(moif/ms)a6, 04, and 04 — 2(mo_fs-/ms)a6, respectively, their real parts decrease in the above order. 
The almost imaginary scalar penguin annihilation then makes the total penguin amplitudes more orthogonal 
to the color-allowed tree amplitudes. As a consequence, the modifications from the NLO color-suppressed 



0.72 
0.67 



0.72 
0.66 



0.71 
0.50 



0.69 0.70 0.84 



+0.01 (+0.01) 
0.01 (-0.00) 
+0.10 (+0.04) 
0.06 (-0.04) 
+0.03 (+0.00) 
-0.07 (-0.02) 
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tree amplitudes decrease in the above order. We then predict the patterns of the direct CP asymmetries in 
Eqs. (fT3|) — (|T5|) . These different patters, if confirmed by the data, will support the source of strong phases from 
the scalar penguin annihilation in PQCD. We emphasize that PQCD is the only theoretical approach available 
in the literature, which can produce the above patterns of the direct CP asymmetries. 

A remark is in order. As mentioned in the Introduction, B meson transition form factors are not factorizable 
in collinear factorization theorem due to the existence of the end-point singularities. Recently, a zero-bin 
regularization for these singularities have been proposed, such that B meson transition form factors become 
factorizable [4(|. The idea is to subtract soft modes from collinear effective fields to avoid double counting of 
soft dynamics. When implemented in a factorization formula, this regularization scheme is equivalent to cutoffs 
at small momentum fractions x. The zero-bin regularization has been extended to annihilation amplitudes in 
two-body nonleptonic B meson decays, which then also become factorizable (4lj . As expected, the annihilation 
amplitudes are real, because internal particles, without carrying parton transverse momenta, do not go on mass 
shell, when one introduces cutoffs at small x. Therefore, whether a physical strong phase is generated by the 
LO B meson annihilation amplitude is a distinction between this modified collinear factorization theorem and 
the kr factorization theorem. 

We have also predicted the mixing-induced CP asymmetries of the penguin-dominated B — ► PV decays. 
Among the decays we have investigated (B° — > n°Ks, <fiKs, p K$, luKs), B° — > (f>Ks is the cleanest mode for 
extracting the standard-model parameter sin(20i), since it does not involve the color-suppressed tree amplitude, 
which can be greatly enhanced by the NLO effects. It has been found that S p o Ks exhibits the maximal negative 
deviation from those of the tree-dominated b — > ccs transitions, and all other Sf show positive deviations in 
NLO PQCD. Hence, it is still puzzling, when confronting the theoretical predictions, especially those for S$k s 
and S 7T o Ks , with the experimental observations. 
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APPENDIX A: DISTRIBUTION AMPLITUDES 



The pion and kaon distribution amplitudes (33j . and the p and K* meson distribution amplitudes [31 ,] have 
been updated up to two-parton twist-3 level in QCD sum rules. The twist- 2 pion (kaon) distribution amplitude 



^n(K) anc ^ the twist-3 ones 4>n(K) an< ^ ^(k) nave been parameterized as [3 



2V2N r 



K(K) r< 3/2 



Cl<*{2x - 1) + al {K) Cl /2 {2x - 1) 



(Al) 



and 



2J2N r 



l + ip 



rr(K) 



27 V (K) ir(K) tt(K) ( 3 



+Cl /2 (2x-l 

+C 1 2 /2 (2x-l) (30r, 3nm + 15p 



tt(K) 



TT(if) TT(if) 

«2 



27a 



+Cl /2 (2x - 1) (l0r] 3n (K)^3n(K) - ^ p^ K) a2 {K) ^J - 3rj 37r (K)U37r(K)Cl /2 (2x - 1) 



7T(K) 



tt(K) 



60. 



tt(K) 



+ | (pl {K) -p^)(l + 3a^+6al 



(j<: 



lax 
lax 



(A2) 
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2J2NZ 



(1 - 2x) 



1 + - p1 {K) 



1( . MK) ir(K) 15 n {K) k(K) 
10 P+ "2 2* 



„ 7T(K) tt(K) 15 7r(lf) n(K) 

°P+ u l F- "2 



1 



+ ^5?/ 37T (K) ~ ^ %7r(iC) w 3ir(K) + g ^+ "2 
+ 5?737r(K) A 37r (/f ) (1 - 2x) [35x(l - x) - 2] 



3 (2a;- 1) 



^[ P \ {K) +P AK) )[1-Za[ 



( 7TI 

2 W 



1 + 3ai 



In a; 

ln(l - ai) 



3p+ (Ar V (in - y ^ (Ar) a2 (if) ^ 6ai(l - x) - l5r] 3n(K) X 3lr{K) x(l - x) 



3 
3 

■2^' 



(p\^+f} K) ) (l-3aj 



>)(!-*) 



7r(iC) _ ^tt(K) 



1 + 3% 



respectively, with aj" = 0, = 0.06, = a% = 0.25 and 

f%n{K) m q(s) + 



■k(K) 
P+ 



P- 



(K) 



rr> 



l(K) U(K) m 07r(K) ' 



(m g(s) + m q ) 2 _ m n{K) 



Ott(K) 



m g( S ) - m \ 



m 



(A3) 

(A4) 
(A5) 
(A6) 



t(K) 



where m g is the mass of the light quarks u and d. One has p+ = and = — for A", and numerically 
pK ^ pK an£ j p!^ ^ q ^j ie ex pii c it values of the other relevant parameters are — J 3 k = 0.0045 GeV 2 , 
A3tt = 0.0, Xzk = 1-6, (^3,1 = —1-5, and uj 3 k = —1-2 [33|. In the above kaon distribution amplitudes the 
momentum fraction x is carried by the s quark (43|. The Gegenbauer polynomials are defined by 



C*(t) = 1 , 

C?(t) = 2At, 

C 2 A (i) = 2A(A + l)i 2 -A, 

-^(A 2 + 3A + 2)t 3 - 2A(A + l)i , 



C£(t) 



C*(i) = ^(A 3 + 6A 2 + llA + 6)i 4 -2A(A 2 +3A + 2)t 2 + -(A+l) 



(A7) 
(A8) 
(A9) 

(A10) 
(All) 



The twist-2 distribution amplitudes for longitudinally polarized vector mesons are parameterized as |3l] 

3/p 



<t>p{x) 
<f><p (x) 



V2Ae 

V^Nc 
Vk- 

x/27^ 



ie(1 


-x) 


1- 


-4 p C 3 2 /2 (2x- 


1) 




x(l 


-x) 


1 - 


~a^Cl /2 {2x~ 


1)" 




x(l 


-x) 


1 - 


^4k* Cl /2 {2x 


-l) + al K »C 3 2 /2 (2x-l) 


a:(l 


-x) 


1 - 


~a^C 3 2 /2 (2x~ 


1) 





(A12) 
(A13) 
(A14) 
(A15) 
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For the twist-3 vector meson distribution amplitudes (f> s and we adopt their asymptotic models [13]. Since 
the Gegenbauer coefficients of these twist-3 distribution amplitudes have not yet been well constrained, the 
asymptotic models are allowed. Moreover, we shall vary the Gegenbauer coefficients of the twist- 2 distribution 
amplitudes by 100%, which is larger than the error specified in [3l|. Therefore, the theoretical uncertainty of 
our predictions from this source is representative enough. 
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